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Abstract 

We study the properties of topological defects in the lattice compact Abelian Higgs 
Model with charge Q = 2 matter field. We find that monopoles and antimonopoles 
form chain-like structures which are dense in the confinement /symmetric phase. In this 
phase the mentioned structures explain both the confinement of single-charged and the 
breaking of strings spanned between doubly-charged test particles. This observation helps 
to understand how the non-diagonal gluons, once taken into consideration in the Abelian 
projection of gluodynamics, could reproduce in this framework the string breaking for 
adjoint charges. 



1 Introduction 

In a series of papers, we have studied the three-dimensional compact Abelian Higgs model 
(cAHM 3 ) as a toy model pQ for confinement and deconfinement in QCD. Starting from three- 
dimensional compact QED (CQED3), where monopoles and antimonopoles in the plasma phase 
are the agents of confinement at all couplings [2j, we have studied deconfinement either in- 
duced by rising temperature (in cQED 2 +i [3) or occurring due to interaction with a single- 
charged (Q — 1) matter field in the framework of the cAHM 3 . In the latter case, the presence 
of matter fields eventually forces monopoles and antimonopoles to form magnetically neutral 
bound states pQ connected by Abrikosov-Nielsen-Olesen (ANO) vortices j3]. This mechanism 
of dipole formation is obviously different from that observed in the high temperature phase of 
cQED 2+1 0E1. 

The role of monopoles is paralleled in four-dimensional Yang-Mills theory where, in the dual 
superconductor scenario jB] , monopoles are also playing the role of agents of confinement (see, 
e.g. the reviews [7]) and where the presence of a dynamical matter field changes the monopole 
dynamics in a way resulting in string breaking -8.j. Not only the similarities, also the differences 
to the four- dimensional non- Abelian theory are interesting in this context. One has to recognize 
that the dynamics of particle-like monopoles in four dimensions is richer and changes at the 
deconfinement transition in a way fundamentally different from the dynamics of instanton-like 
monopoles in three dimensions. 



In 3 + 1 dimensional non-Abelian gauge theory, the Abelian monopoles are condensed at 
low temperatures and electric charges are confined [7]. However, as the temperature rises 
the monopoles become massive and static and the monopole condensate disappears. The 
static monopoles are unable to support confinement of timelike moving charges. However, 
such monopoles can still give rise to a non-zero spacelike string tension. 

In 2 + 1 dimensional Abelian gauge theory, the monopoles form a plasma phase at low 
temperatures in which the electric charges are confined similarly to 3 + 1 dimensional theory. 
Increasing temperature leads to the formation of magnetic dipoles out of monopoles and anti- 
monopoles. At sufficiently high temperatures the monopole plasma is destroyed and a dilute 
dipole plasma is formed. However, the dipole plasma is unable to create neither a "spatial" 
area law nor the confinement of static charges [5J E] • 

In four-dimensional non-Abelian theory there is another, rivalling scenario of confinement 
based on center vortices jU]. A closer look reveals that it does not invalidate the monopole 
picture of confinement but puts it into another context, thus adding essential features to the 
monopole dynamics. In short, percolation or static behavior of the latter results from the 
respective properties of the vortices in the sense of monopoles running inside the vortex sheets. 
In the case of SU(2) this means that monopole worldlines are confined to vortex sheets, and 
magnetic flux between opposite monopoles is collimated inside the vortex sheet. This picture 
has been advocated by Greensite et al. [TU] from the moment that they revived the vortex 
scenario of confinement. Later on, the picture helped to understand the different confinement 
properties (and scales) of fundamental and adjoint charges which were difficult to explain in 
the monopole picture jTTj. 

The cAHM 3 is of interest to mimic this particular aspect due to its non-perturbative features 
which arise from the presence of two types of topological defects: monopoles and Abrikosov- 
Nielsen-Olesen vortices j3]. In the single-charged cAHM 3 tight monopole pairs are bound by 
ANO vortices. In three dimensions monopoles and anti-monopoles experience a logarithmically 
rising attractive potential due to an anomalous dimension of the photon propagator induced 
by the matter fields ^2] ■ The interaction guarantees the formation of monopole-antimonopole 
pairs. When such pairs are forming, asymptotic confinement is simultaneously lost for all test 
charges q. This is not the case we are now interested in. 

In this paper we therefore consider the cAHM 3 with doubly-charged (Q = 2) Higgs fields. 
Besides of its importance in condensed matter physics [13] this model is attractive in our context 
because it is clear [T3j that the dynamical Higgs field can only screen external particles with 
even charge q. Before the transition to the Higgs phase sets in, the q = 1 test charges (the 
analogue of fundamental charges) should be confined at large separations, while q = 2 test 
charges (the analogue of adjoint charges) should suffer string breaking. The assumed formation 
of only magnetic dipole states is not able to explain this selective string breaking. Finally, in the 
Higgs phase, none of the test charges will be confined, which corresponds with the formation 
of monopole and vortex clusters with a size of the order of the lattice spacing. Such clusters 
may consist of pure vortex rings or also they may contain monopoles. 

2 Some analytic considerations 

To enable analytical considerations and to simplify the numerical simulations, we study the 
model in the London limit where the Higgs field $ at site x is represented only by its phase 



$ x = exp (iif x ). The Wilson-type action of the model with Q-charged Higgs field is 

S w [9, <p] = -l3^2 cos(d6) P -kJ2 c os(# - Q0)i , (!) 

p i 

where 9i is the link angle representing the compact gauge field and (dtheta)p is the plaquette 
angle representing its curl. /3 denotes the gauge coupling and k is the hopping parameter. Here 
and below we use the compact notations of differential forms on the lattice (see the second of 
Refs. [7j for a review). 

The monopoles and the ANO-vortices appear due to compactness of the phase angles Q\ 
and ip x , respectively. For the sake of our arguments it is instructive to perform the Berezinsky- 
Kosterliz-Thouless |T5] (BKT) transformation in order to rewrite the partition function of the 
CAHM3 in terms of those defects. We use the Villain-type action (with couplings J3 and k) 
instead of (JTJ): 

Z = [vefVipS^ y^ e -P\\d8+2TTn\\ 2 -K\\d V -Qe+27rl\\ 2 ^ (2) 

-vr _7r ra ( c 2) '( c l) 

where n and I are integer- valued forms living on plaquettes c 2 and links ci, respectively. 

Applying the BKT transformation (T5] (see also the second of Refs. |7[) with respect to the 
gauge and the Higgs fields and integrating them out we rewrite the partition function (J2J) as 

Z oc Z d = J2 J2 z- Sd{ * m * jm) (3) 



*m(*C3) *3m(*C2) 

S*j m =Q*m 



with the defect action 1 



S d (*m, *j m ) = 4vr 2 /3 (*m, — — *m) + Ati 2 ~k (*j m , A + M2 *Jm) • (4) 

~ 1/2 
Here M = Q(k//3) is the tree- level mass of the gauge boson, A is the lattice Laplacian. 

The integer-valued forms *m and *j m represent monopoles and vortices, respectively, living on 

sites and links of the dual lattice. A site *c% (link *C2) of the dual lattice is dual to a cube 

C3 (plaquette c<i) of the original lattice. The constraint in the second sum of Eq. (jSJ) requires 

that the vortices begin (and end) at the monopole (anti-monopole) positions, 5*j m = Q *m. 

This indicates that a vortex carries a fraction 1/Q of the total magnetic flux emanating from 

a monopole. 

Let us now consider the contribution of vortices and monopoles to the potential V q (R) 

between a pair of external test particles with charges ±q, separated by a distance R. The 

potential is given in terms of the average of the Wilson loop W q (R, T) = exp{iq(J,9)} in 

the form V q (R) = — T~ 1 \og(W q (R, T)) + const. The current J runs around a contour of 

rectangular shape, R x T with R <C T. Using the same transformations which led us from 

Eq. © to Eq. (jSJ), the vacuum expectation value of the Wilson loop will emerge factorized 

In the journal version of this article the (A + M 2 ) operator in the first term of the r.h.s. is inadvertently 
written as A -1 . 



as (W q ) = (Wq) h • (W q ) d . The first factor is simply the perturbative self-interaction of the 
external loop 
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due to massive photon exchange. The non-perturbative factor is due to topological defects 
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The first term in the interaction Si nt is a Yukawa-type interaction (which reduces in the limit 
M — > to a usual cQED-like interaction) between a monopole and the charged test particle. 
The second term is given by the linking number IL(*Jq, J) = (*d*jo, A -1 J) 6 ^ between the 
external particle trajectory J and the closed part 2 *jo of the vortex ensemble *j m . The linking 
number describes the Aharonov-Bohm (AB) topological interaction [TEj . 

The contribution of vortices to the potential V q is twofold. At first, and most importantly, 
the closed vortices interact with the electrically charged external particles via the AB effect if 
q/Q ^ /Z . Secondly, the vortices influence the monopole dynamics via the constraint, contribut- 
ing in this way to the potential indirectly, via the monopoles and their interaction with the test 
particles. Since there are Q vortices attached to each (anti-)monopole, the vortices force the 
monopoles to form magnetically neutral states which are dipoles for Q = 1. If Q = 2, however, 
the monopoles may form not only dipoles but also extended loop-like structures depicted in 
Figure H A typical loop structure contains vortex lines which, in general, are non-orientable 







Figure 1: A schematic view of simplest vortex- monopole configurations and a pure ANO vortex 
configuration in the cAHM 3 with Q = 2 dynamical Higgs field. 

closed loops consisting of vortex segments with definite orientations of magnetic fluxes. A junc- 
tion of segments with different flux orientations must contain a monopole. Thus non-orientable 
vortex loops should necessarily contain magnetic monopoles. 

2 A separation of a general vortex ensemble, *j m — *jo + *j' m , into closed vortices, 6* jo — 0, and open ones, 
S*j' m = Q *m, is ambiguous. However, in the sum JSJ) the ambiguity disappears. 



Orientable vortices with magnetic flux cannot lead to confinement of electric charges J7] . 
Indeed, the long-range interaction between such vortices and the Wilson loop may be provided 
only by the AB effect in a form visible in Eq. © (here we implicitly lean upon that we are 
working in a model with a mass gap in which any interactions - except for the topological 
ones like via the linking number - must be suppressed). However, numerical simulations of 
the non-compact version of the AHM 3 show 17] that a confinement phase is absent while the 
AB-interaction does exist. The vortices in this model can only be orientable since monopoles 
are absent. 

On the other hand, in the Maximal Center Gauge of the four-dimensional Yang-Mills theory 
(YM4) the so-called center vortices are known to lead to an area law of the fundamental Wilson 
loop via the linking number representing the AB effect JUJ QB]. The center vortices are non- 
orientable and therefore there is no contradiction with the conclusions of Ref . ^7] . Moreover, 
considered in an Abelian sense the non-orientable center vortices can be represented as segments 
of orientable vortices with monopoles separating the segments with different flux orientations 
from each other, similarly to what we will find in the doubly-charged cAHM 3 . 3 

We conclude that in the doubly-charged cAHM 3 both monopoles and vortices are responsible 
for the confinement of q — 1 charges and the breaking of the string between q = 2 charges. 
Note that in the single-charged CAHM3 non-orientable vortices do not exist, the AB effect is 
absent and the vortices do not feature in producing the confining forces. In this case monopoles 
alone are responsible for confinement while the vortices can only create magnetic dipoles and, 
in this way, destroy confinement. 

3 Numerical simulations 

To check these general ideas, we have performed a numerical study of the confining or non- 
confining properties of electrically charged test particles in the doubly- charged cAHM 3 . The 
phase structure of the model was investigated in Refs. [TUEI]. Before the Higgs phase sets in, 
external particles with charge q = 1 must be linearly confined in this model whereas particles 
with charge q = 2 show a flattening of the potential at a certain distance. In order to observe 
this effect we have simulated the model (0) on a 32 2 x 8 lattice (L t = 8). The choice of the 
asymmetric lattice is dictated by the fact that in the case of symmetric lattices the potential 
can be measured only using Wilson loops which are not suitable for an observation of the string 
breaking effect. We have performed simulations at fixed gauge coupling constant, (3 = 1.2, the 
choice of which was motivated by visualization reasons: to see clearly the monopole structures 
the density of the monopoles must be neither too high nor too low. We have used about 10 4 
independent measurements for each value of the hopping parameter k. 
Defining a Polyakov loop of charge q at position (xi,X2) via 

U 



D q (x 1} x 2 ) = expliq ^ 6 3 (xi,x 2 , x t )\ 



x t =l 



with the obvious notation Q\ = 6 fJi (xi,X2,x t ) we show the v.e.v. of the position average P q in 
Figure I2fa) as a function of the hopping parameter k. At small (large) n the v .e.v 's of both 
loops are low (high) which corresponds to the confinement (Higgs) phase. One can clearly see 

3 Thus, one can ascribe the leading role in the YM4 confinement phenomenon solely to the vortices (as in 
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Figure 2: (a) The q — 1, 2 Polyakov loops and (b) their susceptibilities vs. k. 



that the average Polyakov loop for the q = 1 external charges, which are really confined on 
the left of the transition, is more sensitive to the transition than P^. The much less sensitivity 
of P2 is consistent with the observation of the q = 2 string breaking (due to Q = 2 particles 
popping up in the vacuum) discussed below. 

In Figure I2^b) we show the Polyakov loop susceptibilities vs. k. The susceptibility xp 2 
related to q = 2 external particles changes between roughly constant levels in the two phases 
at the transition point. At the same time, xp± shows a clear maximum signaling the transition 
from confinement to deconfinement. We have fitted the susceptibility xp± i n the vicinity of the 
maximum by the function xp\ = C\j{Ci + (k — n c ) ) with C li2 , 7 and n c being fit parameters 
and we localize the transition at k c = 1.300(1). 

The potentials 



V q (\R\) = --\og(P q (0)P;(R) 



(9) 



extracted from the Polyakov loop correlators are presented in Figure El for both values of the 
electric test charges q — 1,2. We show the potentials both in the confinement (at k = 1.275) 
and the Higgs phase (at k = 1.325). In the confinement phase the potential for the q = 1 
external charges is approximately linearly rising at large distances (the flattening for R — > 16 
is a result of lattice periodicity). The q = 2 potential shows a rapid flattening corresponding to 
the dynamical Q = 2 particle creation from the vacuum and, eventually, to the string breaking 
already in the confinement phase. In the Higgs phase all potentials show flattening due to 
deconfining nature of this phase. 

Let us now focus on the topological defects which should explain this. The simplest char- 
acteristic of a topological defect is its density. The monopole and the vortex densities are 
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Ref. ^3 E|) but i n this case the vortices must be non-orientable which implies the existence of monopoles 
living on the vortex sheets. 
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1,2 potentials in the confinement (k = 1.275) and the Higgs (k 
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respectively. The monopole charge is defined in the standard way, 

m = ^ d [ d #k > 



'ID 



where [■ ■ -] 2 denotes the integer part modulo 2n. The vortex current is, following Ref. [19] . 
defined as 



ZTT 



(12) 



According to Figure0]the densities (shown in lattice units) of the monopoles and vortices are 
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Figure 4: The monopole and vortex densities vs. k. 

gradually decreasing functions of k already in the confinement phase. Towards the Higgs phase, 
the density of monopoles and vortices drops faster, a fact which agrees with the expectation 
that the confining properties of the CAHM3 are due to the topological defects. 



One has to realize that in the confinement phase the monopoles cannot be in a plasma state 
because in this case the potentials for both external charges q = 1, 2 would have to be linearly 
rising at large separations, contrary to our observation (Figure EJ). On the other hand, the 
monopoles cannot form magnetically neutral monopole-antimonopole bound states (dipoles) 
exclusively because in this case both q = 1, 2 potentials V q would have to show flattening - 
essentially for distances larger than the dipole size - again in contradiction to the observed 
behavior (Figure EJ). Thus the only possible kind of monopole configurations - which could 
explain both the linearly rising potential for the q = 1 electric charges and the string breaking 
for the q = 2 charges - is a monopole chain schematically plotted in Figure ^ 

If a vortex loop contains monopoles, the monopoles and antimonopoles are mutually alter- 
nating along the chain. Thus the magnetic flux coming from a monopole inside the chain is 
separating into two parts, gradually squeezing into vortices of finite thickness and, as a conse- 
quence, forming a non-orientable closed magnetic flux. Each piece of such a vortex carries in 
average a half flux, $ = 2-k/Q = ir. If such a flux pierces the g-charged Wilson loop it provides a 
contribution to the loop close to (— I) 9 . This leads to the necessary disorder for odd-charged ex- 
ternal particles (eventually leading to a linearly rising potential) whereas even-charged particles 
are not confined. In Figure El we visualize a typical monopole/vortex configuration observed in 




Figure 5: Example of a monopole/vortex configuration in the confinement phase (at k = 1.275); 
monopoles are shown by circles, vortices by lines. 



our numerical simulations. This example shows the presence of the monopole chains confirming 
the physical picture described above. 

One can suggest that in the confinement phase the monopole chains are percolating as in 
Figure El (so that the monopoles are relevant to infrared physics similarly to the monopoles in 
QCD |7J) whereas in the Higgs phase the monopole chains are relatively short. To check this idea 
we have analyzed the cluster structure and show in Figures Efa,b) the normalized distributions 
of the observed mutually disconnected clusters with respect to the number of monopoles and 
antimonopoles N they contain (Dm) and with respect to the number of vortices links (the 
vortex length L) (Dy). One can see that in the confinement phase, k < K c , the distributions 
Dm(N) and Dy(L) have a two-peak structure. The first peak corresponds to the ultraviolet 
clusters and is centered at small monopole numbers N or small length L, respectively. 4 The 



4 Vortex loops without monopoles obviously are not shown in the histogram of D 
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Figure 6: Histograms of connected vortex clusters (as observed per configuration), (a) with 
respect to the monopole plus antimonopole number N and (b) with respect to the vortex 
length L (the number of links in the cluster), are shown for various values of k. Confinement 
of q = 1 test charges is associated with cluster percolation. 

second peak is due to large (infrared) chain-like clusters which extend all over the volume of 
the lattice. An example from the confining phase near to the transition is shown in Figure EJ 
In the confinement phase the infrared peak in both distributions is clearly present whereas in 
the Higgs phase the infrared peak is absent. This is the signal of deconfinement. 

4 Discussion and outlook 



Summarizing, we have found that in the presence of a doubly-charged Higgs field the monopoles 
must form chain-like structures. This offers an explanation of both confinement of single- 
charged electric particles and string breaking for doubly- charged test particles. This physical 
picture - observed in our study in the compact Abelian gauge theory with a Q = 2 dynamical 
Higgs field - has a close analogy with non- Abelian gluodynamics where tight correlations be- 
tween Abelian monopoles and center vortices (each in the respective Abelian projection) have 
been found [20] • Our observation also suggests a natural way for the formation of monopole 
sheets (instead of chains) in gluodynamics. For example, in the pure SU(2) gauge model (cho- 
sen here for simplicity of discussion) the Abelian monopoles are defined with the help of an 
Abelian gauge, in which the off-diagonal gluons (originally ignored in the Abelian projection) 
play the role of the doubly-charged matter fields coupled minimally to the leading diagonal 
gluons. These matter fields may cause the monopole trajectories to be confined inside sheets, 
a mechanism which in turn should be responsible for the simultaneous occurrence (in the pure 
gauge model) of confinement for fundamental charges (quarks) and flattening of the potential 
between adjoint charges (gluons). A more detailed study of monopole chain formation and 
confining properties of the Q = 2 model will be presented elsewhere [2*Tj . 
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